Hypersound damping in vitreous silica measured by picosecond acoustics 
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The attenuation of longitudinal acoustic phonons up to frequencies nearing 250 GHz is measured in 
vitreous silica with a picosecond optical technique. Taking advantage of interferences on the probe 
beam, difficulties encountered in early pioneering experiments are alleviated. Sound damping at 
250 GHz and room temperature is consistent with relaxation dominated by anharmonic interactions 
with the thermal bath, extending optical Brillouin scattering data. Our result is at variance with 
claims of a recent deep-UV experiment which reported a rapid damping increase beyond 100 GHz. A 
comprehensive picture of the frequency dependence of sound attenuation in «-Si02 can be proposed. 

PACS numbers: 61.43.Fs, 63.50.+X, 78.47.+p 
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Despite considerable recent activity, high frequency 
sound attenuation in glasses remains a controversial is- 
sue. Several processes lead to sound-wave damping in 
glasses. At temperatures T above the quantum regime, 
the most important ones are: thermally activated relax- 
ations [l|, HJ , interaction with thermal vibrations [|| , scat- 
tering from disorder [USB], and resonance with low- 
lying modes [7|, [8|. Sorting these out in function of T 
and the sound-wave angular frequency fl proves to be 
difficult. 

In acoustic experiments, one measures the energy at- 
tenuation coefficient a, or the linewidth of Brillouin scat- 
tering lines r — av, where v is the sound velocity and T 
is in angular frequency. At low Cl and sufficiently low T, 
it is generally agreed that a is dominated by thermally 
activated relaxation (TAR) of structural "defects" . In 
vitreous silica, v-Si02, a strong broad attenuationpeak is 
observed around 50 K in ultrasonic experiments Qj . This 
is well described by the phenomenology of TAR (2j, al- 
though the exact nature of the "defects" is not clear. The 
damping depends on their strength, density, and distri- 
bution. A second damping mechanism, which dominates 
in crystals, is anharmonicity. In the Akhiezer descrip- 
tion, sound modulates the population of thermal bath 
modes which relax via anharmonic interactions 0] . This 
relaxation can be characterized by a mean thermal life 
time Tth- As Tth(T') is generally short, anharmonicity is 
expected to dominate at higher frequencies. It leads then 
to a Brillouin linewidth T an h oc Vl 2 as long as firth << 1- 
This damping contribution increases faster with £1 than 
the Ttar relaxation, as the latter should saturate [l(|. 
For u-Si02, detailed estimates show that at a typical Bril- 
louin light-scattering (BLS) frequency (« 35 GHz) and 
room temperature, TAR and anharmonicity mechanisms 
provide approximately 1/3 and 2/3 of the observed lon- 
gitudinal acoustic (LA) attenuation, respectively [Toj . 

The existence of a third strong damping mechanism is 
evident from the universal plateau in the thermal con- 
ductivity of glasses around 5 K, which implies that at 



high ft the acoustic damping increases dramatically [111 ]. 
Based on the kinetic theory of thermal transport, a must 
increase at least as f2 4 to produce a plateau [12|]. Using 
inelastic x-ray scattering (IXS), direct spectroscopic ev- 
idence for such a dependence, V oc SI 4 , was recently ob- 
tained for the LA waves of two network glasses: densified 
silica [l3| and a lithium borate glass However, IXS 
experiments being in practice limited to scattering vec- 
tors Q > 1 run -1 , the (f2, Q) region where the Brillouin 
linewidth should exhibit this rapid increase remains inac- 
cessible in most cases, including silica. That this regime 
exists also in u-Si02 is fully endorsed by the observation 
that BLS linewidths do not extrapolate smoothly to IXS 
measurements, as shown below. The origin of this atten- 
uation is currently debated. On the one hand, one might 
assume scattering of the sound waves by frozen-in disor- 
der, e.g. by local elastic constant fluctuations On 
the other hand, there can be a resonant interaction of 
sound waves with quasi-local vibrations [8j, [16| . Recent 
quantitative estimates rather support the second mech- 
anism [l7j . In both cases, the explanation links to the 
boson peak, i.e. to an excess of low-frequency modes 
over the Debye expectation, another controversial feature 
of glasses. The above was recently challenged by deep- 
UV Brillouin scattering experiments claiming an unex- 
pected sharp increase in the acoustic attenuation of sil- 
ica at room-T already around 100 GHz [TU, i.e. much 
below the boson peak frequency of ~ 1 THz. These un- 
settled questions clearly call for more sound-absorption 
data between 100 GHz and 1 THz. 

In this Rapid Communication, we report new measure- 
ments of the LA attenuation in v-Si02 thin films in this 
crucial frequency region using a sensitive modification 
of the picosecond ultrasonic technique (PU) pioneered in 
the 1980s Difficulties encountered in the early work 
(2fl [2l| can be alleviated using the present data acqui- 
sition and analysis. Specifically, we find damping coeffi- 
cients approximately one third those in [20| , in line with 
available BLS data. The relaxation via anharmonic in- 
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FIG. 1: (a) Schematic diagram of the experiment; (b) Re- 
flected intensity as a function of time delay obtained on an 
Al/SiO2(600 nm)/Si sample using a 820 nm pump. 



teractions with the thermal bath appears to be the dom- 
inant broadening mechanism at 300 K for £1/2it between 
- 190 and ~ 240 GHz. 

The samples are v-SiC>2 layers of four thicknesses from 
300 to 1200 nm deposited on a (100) silicon substrate 
by low pressure chemical vapor deposition. This pro- 
duces layers with macroscopic properties very similar to 
those of bulk fused silica. The layers contain ~ 3000 
ppm OH plus a small amount of hydrogen, below 500 
ppm. An Al film of thickness d ~ 12 nm is evaporated 
on top of the samples in a single batch, to obtain iden- 
tical films on all of them. We use the standard pump- 
probe setup sketched in Fig. 1(a). The light source is a 
titanium:sapphire oscillator tunable in the 700-980 nm 
wavelength range, producing 120 fs pulses. The pump 
is directly derived from the oscillator, whereas the probe 
is frequency doubled and time delayed. The pump and 
probe fall on the sample surface at near normal incidence, 
focused to a spot of ~ 20 fim diameter. The pump is 
absorbed by the thin Al film which expands launching 
a short longitudinal strain pulse into the sample. The 
generated acoustic pulse is a plane wave packet of ~ 20 
nm spatial extent whose Fourier spectrum peaks around 
250 GHz, following 22]. This pulse undergoes a reflec- 
tion at the layer/substrate interface producing an echo in 
the Al film which can be detected by the delayed probe. 
In the early work by Zhu, Maris, and Tauc [20j| . the fre- 
quency dependence of the Si02 layer attenuation was de- 
termined from the decay of successive echoes. A tung- 
sten substrate was used to produce numerous and str ong 
echoes. Here, we use an interferometric method [23], [24[ 
to compare the strength of high frequency components of 
the acoustic pulse reaching the Si substrate in function 
of the layer thickness [H, ^ H3 . 

The acoustic pulse is detected by means of the time- 
delayed probe as indicated in Fig. 1(a). The main part 
of the probe is reflected at the fixed interfaces. A weaker 
component is reflected by the moving strain pulse, ac- 
cording to the Bragg relation, q = 2k cos 9 [23| |. where q 
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FIG. 2: (Color online) (a) Typical fits (solid lines) of the sig- 
nal intensity; (b) Intensities normalized to the A\ amplitude, 
labeled by the nominal Si02-layer thickness; (c) The normal- 
ized ^b-value as a function of the Si02-layer thickness. 



is the wave vector of the LA phonon, and k and 9 are the 
wave vector and the angle of incidence of the probe light, 
respectively, both inside the material. The two reflec- 
tions interfere leading to oscillations with the time delay 
t. The period is 

T = A probc / (2nv cos 9) , 

where A pro b e is the probe wavelength, n is the refractive 
index at A pro be, and v is the LA sound velocity. 

A typical signal is shown in Fig. 1(b). At zero delay, 
the intensity exhibits a sudden step followed by a long 
decay related to the thermal response of the Al film. Up 
to the transit time t\ = 98 ps, a sinusoidal modulation 
of long period, T\ — 23.5 ± 0.1 ps is superposed to this 
background. With n = 1.470 [HI, T x leads to the correct 
LA sound velocity of i?-Si02, v = 5940 ±25 m/s. Beyond 
t%, one observes strongly damped oscillations of smaller 
period, T-i = 4.57±0.02 ps, in accordance with the refrac- 
tive index 5.32 [29| and the LA sound velocity 8445 m/s 
[30l | of the Si(100) substrate. The frequency of the Fourier 
component probed in the Si02 layer is 1/Ti = 42.5 GHz 
while it is much higher, 1/Tq, = 219 GHz, in the silicon 
substrate. One also observes during the first picoseconds 
a superimposed damped oscillation. This is the strain 
pulse partially bouncing back and forth within the Al 
film. From its transit time, 1.9 ± 0.1 ps, we calculate 
an Al-film thickness of 12.2 ± 0.6 nm, the same for all 
samples. This is an important point as the frequency 
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FIG. 3: (Color online) Full widths T(f2) of the Brillouin peaks 
deduced from different spectroscopies, as indicated. The full 
line is the estimated sum of TAR (dotted line) , anharmonic- 
ity (dashed line) and SPM (dot-dashed line) contributions as 
explained in the text. The thin solid line of slope 2 is a guide 
for the eye. The inset illustrates the unexpected sudden rise 
of linewidth at ~ 100 GHz reported in Ref. 18. 



distribution of the acoustic pulse is controlled by this 
thickness and our analysis assumes identical pulses in all 
samples. 

For the analysis the background is first subtracted from 
the raw data. Typical signals after subtraction are shown 
in Fig. 2(a,b) [31]. The oscillations appear undamped 
in W-S1O2 because both the light and sound absorptions 
are negligible at those frequencies. In contrast, a strong 
damping rate of the oscillations, 5.51 x 10~ 2 ps~ 4 , is 
observed in silicon, in agreement with the high light ab- 
sorption coefficient of Si at A pro be [HI- We call A\ and 
A2 the initial oscillation amplitudes in the SiC>2 layer 
and in the Si substrate, respectively. The acoustic at- 
tenuation of the layer at the high frequency, I/T2, can 
be obtained measuring the decrease of A2 with increas- 
ing layer thickness, as illustrated in Fig. 2(b)-(c). For 
normalization, the signal divided by Ai is plotted in 
Fig. 2(b). To correctly normalize A2 one must include in 
Ai both the Stokes and anti-Stokes processes, plus the 
effect of the Fabry-Perot interferometer formed by the 
t>-Si02 layer which affects both Ai and A2. Finally, the 
normalized A2 value follows an exponential decay with 
increasing Si02 thickness with an attenuation constant 
1.80±0.15 /im -1 , from which we deduce the energy mean 
free path, t = a -1 , of LA phonons in the silica layer to 
be 280 nm at 219 GHz and room temperature. 

Tuning the laser wavelength we explored the LA 
phonon attenuation over the range 190-250 GHz. Our 
results are shown in Fig. 3 together with the early PU 
data [20ll . and with the available BLS [H, HU, [34j and 
IXS [35j| data. To compare these results all are presented 
as Brillouin full widths. Also shown in Fig. 3 are the 
predictions of the most widely accepted models for the 
three main damping mechanisms. As explained above, 



the thermal activation of relaxing defects (TAR) accounts 
for the sound attenuation at sonic and ultrasonic frequen- 
cies while the network viscosity or anharmonicity takes 
over at the BLS frequencies. The dotted and dashed lines 
in Fig. 3 illustrate the expected contribution of the TAR 
and the anharmonicity broadening mechanisms, respec- 
tively, following [10], on the basis of experiments covering 
7 decades in frequency [36| . Here, the curves are extrap- 
olated up to the THz frequency range. Moreover, one 
expects in this frequency range the onset of the sound 
attenuation oc f2 4 discussed above. This strong growth 
can lead rapidly to the end of acoustic plane waves at the 
Ioffe-Regel limit, which does provide an explanation for 
the universal low-T plateau in the thermal conductivity 
pd| . In u-Si02 this limit is anticipated at Qjb ~ 1 THz, 



just at the boson peak maximum 37J. It corresponds to 
a wave vector q ~ 1 nm -1 for the LA mode, a region 
which is not attainable with current IXS spectrometers. 
The full square in Fig. 3 obtained at 1.5 nm -1 is, in fact, 
the lowest accessible IXS point measured at room-T in 
t>-Si02 [35j]. Within the framework of the soft poten- 
tial model (SPM) (38[, a quantitative prediction for the 
L oc fl 4 behavior can be calculated [3£|, shown as the 
dash-dotted line. This f2 4 dependence is expected to ter- 
minate at the Ioffe-Regel limit, indicated as an arrow in 
Fig. 3. The solid line is simply the sum of the three 
separate contributions. 

In the low frequency part, it can be seen that the dif- 
ferent fl and T dependences of Ttar and r an h lead to 
a non-trivial variation of L with Q. The departure from 
the f2 2 behavior is indeed evident for the BLS results over 
the frequency range 20-100 GHz. At higher frequencies, 
our PU data become dominated by the f2 2 dependence 
of r an h, in remarkable agreement with the prediction. 
The early PU data follow roughly the predicted depen- 
dence shown by the solid line, although it is obvious now 
that the absolute value was systematically overestimated. 
Several sources of error are in fact discussed in [2(| , the 
most serious one being possibly the corrections for losses 
at the interfaces. The above scenario was challenged in a 
recent Letter reporting deep-UV Brillouin scattering data 
in which a dramatic increase of linewidth is observed in 
the narrow frequency range 100-130 GHz as shown in 
the inset of Fig. 3 The authors interpreted it as the 
onset of a new acoustic damping at a scattering vector 
Qi ~ 0.11 nm -1 up to Q2 — 0.15 ran -1 . They associ- 
ated it with Rayleigh scattering by disorder at the huge 
length scale, 27r/Q 2 - This is surprising as it is gener- 
ally thought that the scale of disorder in a glass does not 
exceed a few nanometers. It is now clear from our re- 
sults that no such damping onset occurs up to at least 
250 GHz at room temperature in i;-Si02. The increase 
of the spectral linewidths observed in the deep-UV Bril- 
louin scattering measurement could somehow be related 
with the optical absorption edge of w-Si02, in which case 
it is not an acoustic attenuation. 

Above the Ioffe-Regel limit, there also exist IXS data, 
e.g. in 35]. As shown in Fig. 3, anharmonicity damping 
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saturates in that region so that it cannot explain the val- 
ues observed in IXS giving an indirect evidence for the 
rapid growth of attenuation. The onset of the T oc il 4 
should occur around 400 GHz at room temperature in 
silica. Above £lm, the exact nature of the excitations, 
which follow an apparent damping r oc f2 2 , is so far un- 
known. The fact that IXS spectra are observed and can 
be approximately adjusted to a damped harmonic oscil- 
lator lineshape is definitely not a proof that these ex- 
citations are propagating plane waves. As discussed in 
pol ]. elastic heterogeneities are seen in silica simulations 
at the boson peak scale. A likely view is that beyond the 
Ioffe-Regel limit the acoustic modes strongly mix with 
quasi-local vibrations and thereby become diffusive, 
as supported by other simulations [411 ]. 

In conclusion, we report new sound-attenuation coef- 
ficients in vitreous silica in the crucial frequency range 
between 100 GHz and 1 THz. These new data points 



prolongate smoothly lower frequency Brillouin light scat- 
tering information. The results fall on a line calculated 
independently on the basis of the three main damping 
contributions active in this frequency-temperature do- 
main. We observe that beyond 100 GHz up to at least 
250 GHz it is anharmonicity that dominates sound damp- 
ing at 300 K in silica. This implies that the regime of 
strong growth of the attenuation should arise beyond this 
frequency range, exactly on the low frequency side of the 
boson peak in silica. Thanks to the present experimental 
breakthrough one can now envisage to observe the onset 
of the expected f2 4 dependence. As the latter mechanism 
should be T-independent, the onset is expected to move 
down to an accessible frequency range as T is decreased 
reducing the anharmonicity damping. This will be a fu- 
ture crucial test of the overall picture. 
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